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ABSTRACT 


The phenomenon of the rapid growth and development of an 
extratropical cyclone over the east coast of the United 
States (the Carolinas storm of March 1984) is studied 
through a linear stability analysis. Analyses of the 
cyclone structure suggest barotropic and baroclinic insta- 
bilities may be important. A linear stability model is used 
to investigate the roles and relative importance of short- 
wave baroclinic instability and barotropic instability in 
the growth and development of the storm. The growth rates, 
phase speeds and structure obtained from the linear model 
are consistent with those derived from observations. Energy 
budget results indicate that the vertical and horizontal 
barotropic terms are at least as important as the baroclinic 
term. It appears that the early growth and development of 
the cyclone can be explained through the contributions of 
barotropic and baroclinic instabilities without including 


convection. 
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I. INTRODUCTION 


Extratropical cyclones have been the focus of extensive 
research for many years. Modern studies of genesis and 
development of cyclones began with the pioneering work of 
the Norwegian school that formulated the polar front theory 
(Bjerknes, 1919). Later, Charney (1947) proposed the insta- 
bility theory to explain the cyclogenesis process. 
Petterssen (1956), Palmen and Newton (1969), Holton (1979) 
and others have provided textbook treatments of the struc- 
ture and behavior of a typical cyclone. 

A number of factors indicate that the east coast of the 
United States is a primary location for Significant cyclone 
development (e.g., Sanders and Gyakum, 1980). These include 
land-sea temperature contrasts, the location and intensity 
of the Gulf Stream, air mass modification over the Atlantic 
Ocean, land-sea frictional effects, the shape of the coast- 
line and the influence of the Appalachians (Kocin and 
Uccellini, 1985). Some of the cyclones that form along the 
east coast deepen rapidly and have presented major forecast 
problems. Operational as well as research numerical fore- 


casts have only had marginal success in predicting the deep- 


ening process (Sanders and Gyakum, 1980; Bosart, OS 
Atlas, 1984). These cyclones exhibit deepening rates of at 
least 1 mb/h for 24 h in sea-level pressure. Observational 


studies such as that of Sanders and Gyakum (1980) have docu- 
mented the climatological aspects of these explosively deep- 
ening cyclones. Deepening rates of these cyclones cannot be 
explained from quasi-geostrophy alone (Sanders, 1971; 
Sanders and Gyakum, 1980). Observations also have suggested 
that some of them initially form on the meso-alpha scale 
(200 - 2000 km). Therefore, it appears that other physical 


processes may contribute significantly to the growth and 


development of this class of extratropical cyclones. It is 
important that these storms be studied and understood since 
public perception of forecast accuracy is particularly acute 
during these events. These cyclones have a large impact on 
centers of population and business as well as on naval 
operations. 

An example of a rapidly developing cyclone along the 
east coast is the Carolinas storm of 28-30 March 1984. This 
cyclone formed along the North/ South Carolina border with a 
diameter of approximately 1000 km. From 0000GMT through 
1500GMT on 29 March, the cyclone rapidly deepened while 
growing to a diameter of approximately 4000 km. The cyclone 
reached its lowest central pressure of 965 mb at 1500GMT 
while moving northeastward away from the coast. Its central 
pressure fell 15 mb in this 15 h period. 

The operational numerical models (the Limited Fine Mesh 
Model (LFM) and the Navy Operational Global Atmospheric 
Prediction System (NOGAPS) ) failed to predict the rapid 
development of this cyclone. However, the Navy Operational 
Regional Atmospheric Prediction System (NORAPS) and the 
Mesoscale Atmospheric Simulation System (MASS) did forecast 
the formation of the cyclone, although in both models the 
predicted low pressure center is southwestward of that 
analyzed by the National Meteorological Center (NMC). 

The inaccurate forecasts of this and other rapidly 
developing cyclones suggest that certain dynamical features 


Or processes associated with these cyclones may not be prop- 


erly simulated or diagnosed. Bosart (1981), Uccellini et 
al. (1981, 1984) and Sanders and Gyakum (1980) have noted 
certain model deficiencies, including inadequate boundary 


layer and cumulus parameterizations, coarse vertical resolu- 
tion and poor jet streak simulations. However, recent simu- 
lations of these storms using mesoscale models have produced 


Some encouraging results (Kocin et al., 1984). 


The objective of this thesis is to better understand the 
genesis and rapid development of maritime extratropical 
cyclones by studying the case of the Carolinas storm. 
Although over land, this storm exhibited many typical char- 
acteristics of a rapidly intensifying maritime cyclone. The 
abundance of land observations allows the synoptic as well 
as the sub-synoptic scale features to be better understood. 
A linear stability analysis will be performed to determine 
the roles and relative importance of shortwave baroclinic 
instability (Staley and Gall, 1977) and barotropic insta- 
bility in the genesis and subsequent rapid development of 
this cyclone. 

A literature survey will be presented first on the roles 
of baroclinic instability, barotropic instability, jet 
streak dynamics and convection on rapid cyclone development. 
Next, a detailed synopsis of the synoptic events that char- 
acterized the Carolinas storm will be summarized. The 
linear stability model developed by Gall (1976c) will be 
used to investigate the role of baroclinic and barotropic 
instabilities in the genesis of this storm. Finally, 
conclusions and possible avenues for future research will be 


suggested. 


II. REVIEW OF THE MECHANISMS FOR GENESIS AND DEVELOPMENT OF 
E MRS SIRATROBICAL CYCLONES = 


A. BAROCLINIC INSTABILITY 

It is generally agreed that extratropical cyclones on 
the synoptic scale derive most of their energy from the 
release of the available potential energy inherent in air- 
mass contrasts. These disturbances depend on a pre-existing 
baroclinic structure of the basic current. However, the 
kinetic energy of the disturbance can also be obtained 
through the conversion of the kinetic energy of the basic 
current. 

Baroclinic instability (BCI) is a mechanism that has 
been proposed to explain how potential energy of the basic 
cn enbican be converted to kinetic energy of the distur- 
bance. The principles of BCI were introduced by Charney 
(1947) and Eady (1949). Their Studies show that the most 
baroclinically unstable waves have a wavelength of approxi- 
mately 4000 km, which corresponds to that observed to have 
the largest amplitude in the extratropical cyclone at mid 
and upper levels (Holton, 1979). 

Historically, the dynamics and thermodynamics of a 
developing baroclinic disturbance have been explained using 
quasi-geostrophic theory. This theory is based on mid- 
latitude, synoptic-scaled disturbances. Furthermore, the 
atmosphere is assumed to be frictionless and adiabatic. 
Also, the atmosphere can only undergo changes in such a way 
that geostrophic and hydrostatic balances are maintained. 
The theory greatly simplifies the system of equations that 
govern the dynamics and thermodynamics of the atmosphere, 
while preserving the fundamental physics which the original 


equations represent. 


Extending the concepts of this theory, Petterssen 
(1956), Palmen and Newton (1969) and Holton (1979) have 
shown that when all other factors are held constant, the 
process of cyclone development can be accelerated if the 
static stability is reduced and/or diabatic heating (latent 
heat release, air-sea fluxes, radiation) is introduced. 
These other factors include warm air advection and the 
advection of positive vorticity at mid-tropospheric levels. 
As a result, the "self-limiting" process, in which the 
effect of temperature advection is partially offset by the 
adiabatic temperature changes caused by vertical motion, is 
reduced. 


Other observational studies (Sanders and Gyakum, 1980) 


have documented the explosive development of certain 
cyclones. These incipient disturbances are relatively 
shallow and confined mostly to lower levels (Bosart, 1981). 


Several studies have attempted to explain this structure by 
modifying the early baroclinic theories which were based on 
quasi-geostrophy. Mansfield (1974) modified the Eady model 
to include a shallow baroclinic zone in the lowest 1.6 km. 
The results indicate that reducing the static stability in 
the lowest layers shortens the wavelength of maximum insta- 
bility to 1000 km. Duncan (1977) and Staley and Gall (1977) 


found similar results. 


Using a linearized general circulation model, Gall 
(1976a) demonstrated that disturbances with a zonal 
wavenumber of 15 exhibit the maximum growth rate. However, 


results from the nonlinear version of the model suggest 


waves with wavenumbers 5-7 to be the most unstable. 
Therefore, although shortwaves are the most baroclinically 
unstable, nonlinearity modifies the spectrum such that 
longer waves eventually dominate. Inclusion of wave-zonal 
flow interaction (Gall, 1976b) increases the low-level 


Static stability because of the northward and upward eddy 


10 


transport of heat. This retards the growth rate of the 


shorter waves since — they are primarily low-level 
disturbances. Gall (1976c) showed that although waves 
of wavenumber 15 are mostly confined near the surface, the 


vertical extent could be increased by including latent heat. 
In another numerical study, Gall et al. (1979) found the 
most unstable waves to be initially at wavenumber 13. 
However, friction, an increase in static stability and 
destruction of the north-south temperature gradient weakened 
these high wavenumber disturbances with time, so that 
wavenumber / became the most energetic. 

These numerical studies as well as earlier case studies 
of extratropical cyclones (Bjerknes, from Godske et al., 
1957; Palmen and Newton, 1969) have focussed on the evolu- 


tion and the attendant structural changes of the cyclones on 


the synoptic scale. These findings suggést that BCI is a 
dominant mechanism in the genesis of a cyclone. However, 
once the cyclone has developed to the mature stage, the 


vertical tilt of the low pressure trough that was prevalent 
at the genesis stage becomes less apparent. This suggests 
that other processes besides BCI may have increased in 
importance. : The next section will review studies which 
relate both baroclinic and barotropic instabilities to the 


genesis and development process. 


B. BAROTROPIC INSTABILITY . 
As previously discussed, much of the current theory of 
eddies in the atmosphere is based upon the dynamic insta- 
bility of zonal currents (Charney, 1947; Eady, 1949 and 
Kuo, 1949). These studies have shown that zonal flows with 
vertical or horizontal shears typical of the atmosphere were 
unstable to small perturbations. For mathematical reasons, 
these earlier studies involved essentially two-dimensional 
atmospheres in which the wind shears were either entirely 


baroclinic or barotropic. The previous section discussed 


PE 


those studies that focussed on baroclinic instability. This 
section will address the effect of barotropic instability on 
the cyclogenetic process. 

Textbook treatments of barotropic instability are 
contained in Holton (1979) and Pedlosky (1982). They both 
showed that the necessary condition for a region to be baro- 
tropically unstable is that the gradient of the absolute 
vorticity of the mean current must vanish somewhere in the 
region. Furthermore, a sufficiently large value of the beta 
term (which represents the change in the Coriolis force with 
respect to latitude) can barotropically stabilize the 
current. Pedlosky (1982) suggested that baroclinic insta- 
bility is the favored mode for horizontally broad zonal 
currents, whereas very narrow currents would tend to be 
barotropically unstable. These authors further showed that 
barotropic instability converts kinetic energy from the 
basic current to the perturbation, in contrast to the energy 
flow description associated with baroclinic instability. 

Studies using more realistic atmospheric flows with both 


horizontal and vertical shears were undertaken by Brown 


(1969) and Sone Aa Their results showed that baro- 
clinic waves are barotropically damped. At longer 
wavelengths, they found that barotropic waves could exist. 


A later study by Gent (1974) found that waves with the 
largest meridional wavelength are the most barotropically 
unstable. 

Most of these studies constructed mathematically simple 
wind fields, where the horizontal shear was invariant with 
height. A more realistic representation of a wind field was 
used by Gall (1976a). However, using a 40 m/s jet, he found 
that waves of all wavelengths receive most of their energy 
through baroclinic processes. For waves’ shorter than 
wavenumber 15, baroclinic instability is the most important 


contributor to the development of the perturbation. On the 
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other hand, the barotropic terms for the most part hamper 
wave growth. 

Pedlosky (1982) suggested that in zonal currents where 
both baroclinic and barotropic instabilities are possible, 
the energy-transfer characteristics of the most unstable 
wave depend upon the detailed distribution of the zonal 
velocity and the potential vorticity of the basic state. 
Changing these parameters of the jet can profoundly affect 
whether baroclinic or barotropic mechanisms are favored for 
instability, as well as the wavelength of maximum 
instability. 

These ideas were incorporated into a numerical study on 
polar lows conducted by Hodur (1984). He found the presence 
of both barotropic and baroclinic instabilities during the 
genesis of the polar low. Inņnsparticulars by varying the 
intensity of the jet from 40 m/s to 60 m/s, he showed that 
for higher wavenumbers and a stronger jet, the perturbation 
can extract as much energy from eddy available potential 
energy as from zonal kinetic energy. For the weaker jet, 
the baroclinic term became the dominant contributor, which 
is consistent with the simulations of Gall (1976a). In 
Chapter 5, these ideas of how energy for the perturbation in 
the Carolinas storm is extracted from the mean flow through 
baroclinic and/or barotropic processes will be investigated 
further. 

There have not been many studies on cyclone development 
through the mechanism of barotropic instability. This 
thesis will be the first to study this mechanism in conjunc- 


tion with baroclinic instability in a real-case cyclone. 


C. JET STREAK DYNAMICS 

Palmen and Newton (1969) showed that appreciable upper- 
level divergence and vertical motion fields characterize 
synoptic-scaled disturbances in regions of the jet stream. 


Families of synoptic disturbances are typically found along 
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the polar front jet stream. From the study of many cases of 
cyclone life cycles, Petterssen (1956) formulated the rule 
that cyclone development occurs when and where a region of 
positive vorticity advection (PVA) at the upper levels 
becomes superimposed over a low-level baroclinic zone. This 
region of PVA is typically associated with the jet stream. 
Other studies have attempted to test the hypothesis that 
when the static stability is weak, the superposition of 
frontal systems in the upper levels on the planetary 
boundary layer can result in favorable conditions necessary 
for rapid cyclone development. Palmen and Newton (1969), 
Kocin and Uccellini (1985) and others have suggested that 
jet streaks imbedded within upper-level trough/ridge systems 


enhance and focus upper-level divergence, energy exchange 
and momentum transport, all of which can play an important 
role in cyclone genesis and development. Shapiro (1982) 


suggested a way in which upper-level and lower-level fronts 
can be coupled via ageostrophic circulations associated with 
each feature. A pattern of intersecting upper and lower 
level jets has been shown by Uccellini and Johnson (1979) to 
be a classic configuration for triggering the release of 
convective and potential instabilities. The resultant 
latent heat release could then enhance the growth and devel- 
opment of the cyclone. 

A case study by Calland (1983) of rapid cyclone develop- 
ment that occurred southeast of Japan in January, L19793 
showed the importance of the jet streak. Specifically, his 
results suggested that upper-level forcing may play a 
greater role in the initiation of explosive development than 
was proposed by Petterssen (1956). They also point to the 
importance of boundary layer and convective processes during 
the explosive development phase. 

From budget statistics, Conant (1982) illustrated the 


importance of the polar jet streak and an amplifying 
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mid-level trough when the Presidents’ Day cyclone explo- 
Sively deepened. In a study of a North Atlantic polar low, 
Wash and Cook (1985) found that the average horizontal 
advection of vorticity dominated the forcing terms and coin- 
cided with the intensification of the cyclone mass circula- 
tion, the large production of vorticity by the low-level 
convergence and the largest decrease in the sea-level 
pressure. 

Uccellini et al. (1984) identified three jets that were 
deemed important in the development of the Presidents’ Day 
cyclone. Diagnostic calculations revealed an increasingly 
unbalanced flow in the subtropical jet (STJ) associated with 
a super-geostrophic flow. The upper-level divergence asso- 
ciated with the STJ was linked to the intensification of the 
low-level jet along the east coast. Subsequent studies such 
as that by Uccellini et al. (1985) focussed on the ampli- 
fying polar jet. This latter study illustrated the role of 
dynamically-forced ageostrophic vertical circulations on the 
deformation of the tropopause. The dramatic tropopause fold 
which occurred during the trough amplification was consis- 
tent with the subsidence expected from geostrophic deforma- 
tion patterns associated with the jet streak. 

Using an analytic formulation, Hoskins and Bretherton 
(1972) were able to cause frontal formation near the tropo- 
pause and the associated tropopause folding through geos- 
trophic confluence. It is apparent from arguments of 
Shapiro (1982) that a correct combination of stretching and 
Shearing deformations can cause tropospheric folding. The 
resultant intrusion of stratospheric air, as characterized 
by high potential vorticity into the troposphere, could then 
spin up the lower levels and cause genesis. It seems reason- 
able that upper-level forcing would more strongly affect 
changes in the lower layers if the intermediate layers were 


convectively unstable. In this way, convective overturning 
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would efficiently transfer upper-level momentum to the 
surface. This apparent contribution to cyclogenesis through 
convective processes will be discussed further in the next 


section. 


D. CONVECTION 
It was observed by Petterssen (1956) that cases of 


moderate and heavy convection over a large area probably 


will not occur, even if surface heating is intense, unless 
the convection is associated with a cyclone. Palmen and 
Newton (1969) summarized from observations that although 


convective clouds cover very large regions of synoptic 
disturbances, areas occupied by clouds of considerable 
vertical development are small. However, Since most of the 
rising motion associated with the cyclone takes place in the 
convective system, Palmen and Newton (1969) suggested that 
convection and the cyclone are dynamically linked. They 
further argued that the static stability, which measures the 
atmospheric susceptibility to convection, typically acts 
against cyclone development by lessening the effect of 
thermal advection. However, this damping term is consider- 
ably smaller for a given vertical motion when extensive 
condensation within the cyclone exists. If this condensa- 
tion is convective in nature, development can proceed 
further and the cyclone can achieve a greater intensity. 

The generation of available potential energy has been 
shown to be enhanced when precipitation occurs in the warm 
sector of a cyclone (Bullock and Johnson, 1971). Tracton 
(1973) stressed the importance of convection in enhancing 
cyclone development during the developing stage. Sanders 
and Gyakum (1980) noted that strong low-level baroclinity, a 
mobile 500 mb trough and intense convection are associated 
with rapid maritime cyclone development. 

In a study of the Queen Elizabeth II (QEII) storm, 
Gyakum (1983a) found that the storm originated as a shallow 
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baroclinic disturbance, and then developed in response to 
low-level warm advection. He hypothesized that the 
mesoscale conditions during the incipient stages of the 
storm were linked to its subsequent explosive development on 
the synoptic scale through a significant increase in poten- 
tial vorticity. This increase was proposed to be related to 
the deep convection, which appeared to be important in the 
explosive deepening phase. The storm formed on the anticy- 
clonic side of the jet, suggesting that inertial instability 
might also have contributed to the development of the 
cyclone. 

The second part of his study (Gyakum ,1983b) examined 
the dynamic and thermodynamic structure of the QEII storm. 
In this case, BCI appeared to be directly responsible for 
only a small part of the cyclone development. BCI was 
important in that it organized convective effects in such a 
way as to produce a positive feedback between convection and 
BCL. Gyakum was able to establish both observationally and 
numerically the critical importance of cumulus-induced 
heating effects in determining whether a cyclone explosively 
develops. However, key questions remain as to how the 
organization of the cumulus’ clouds are affected by the 
cyclone as well as the important aspects of the interaction 
between the cumulus clouds and the larger-scale features. 

Liou and Elsberry (1985) from heat budget studies over 
the western North Pacific showed that increasing the heating 
rate (sensible and latent) reduces the sea-level pressure of 
a cyclone. This suggests that although diabatic processes 
might weaken horizontal temperature gradients and frontal 
zones, warming of the air column near the area of genesis 


aids in the process of cyclone formation. 


E. OTHER MECHANISMS 
Besides the four mechanisms described above, other 


possible hypotheses have been proposed to explain the 
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genesis and development of extratropical cyclones. For 
example, Petterssen (1956) showed the importance of orog- 
raphy by identifying distinct frequency maxima of episodes 
in both winter and summer of genesis to the east of the 
highest parts of the Rocky Mountains and the Alps. Palmen 
and Newton (1969) suggested that the overall process of 
genesis to the lee of mountains is one in which vorticity is 
first generated ina low-level trough, held fixed to the 
eastern slope and finally overtaken by the divergence region 
associated with an upper-level trough. The rapidity of 
development is at least partially due to the availability 
of significant amounts of vorticity at low levels. 

Recent studies have also pointed to the importance of 
certain boundary layer processes in the genesis and develop- 
ment processes. Analyses of surface observations by Bosart 
(1981) of the Presidents' Day snowstorm suggested thát the 
incipient cyclone originated as a shallow, baroclinic 
boundary layer disturbance along a coastal front in the 
Carolinas. Differential heating, moistening and roughness 
between land and sea in conjunction with cold air damming to 
the east of the Appalachians were proposed to be important 
physical mechanisms that contributed to the cyclone devel- 
opment. Rapid deepening of the cyclone occurred when it 
came under a favorable southwesterly flow ahead of a promi- 
nent, upper-level trough which had moved eastward. from the 
Ohio valley. Energy, moisture and vorticity budgets 
analyzed by Bosart and Lin (1984) showed the importance of 
planetary boundary layer processes and cold air damming. 
However, a recent numerical study (Guo and Hoke, 1985) 
suggests that the latent heat release occurred in the 
Presidents' Day cyclone was more important than the sensible 


heat flux from the North Atlantic Ocean. 
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E. DISCUSSION 

Many ideas resulting from various types of studies have 
been proposed to explain the behavior of a certain class of 
extratropical cyclones. These studies have primarily 
focussed on the mechanism of baroclinic instability, with 
the viewpoint that other processes can augment this mecha- 
nism and thereby explain the behavior of those cyclones that 
rapidly deepen. 

Specifically, baroclinic instability has been used to 
explain the structural and developmental characteristics of 
the typical synoptic-scaled cyclone. However, studies by 
Gall (1976a) and Hodur (1984) suggested that barotropic 
instability also could be an important contributor to the 
rapid development of some of those cyclones. In addition, 
case studies have pointed to other processes that might be 
contributing to the genesis and developmental processes of 
these cyclones. Ihese processes minelude (but are not 
limited to) convection, boundary layer processes and jet 
streak dynamics. 

Each of these additional processes explain to a certain 
extent the behavior of a cyclone that rapidly develops. 
They also provide some physical insight as to why certain 
cyclones behave differently. The fact that the incipient 
cyclone can exist on the meso-alpha scale suggests that 
interaction with the mesoscale might be important. These 
ideas raise further questions about these cyclones. For 
example, why do some of these processes become important for 


a cyclone that rapidly develops as opposed to one that does 


not? Do these other processes also contribute to the 
genesis of the cyclone? If so, how does their combined 
effect modulate the developmental process? Does a strongly 


baroclinic atmosphere provide the background that is condu- 
cive to enhancing the effect of these other processes, or 
vice versa? What additional processes might be important at 


the genesis phase? 
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Sanders and Gyakum (1980) have shown that episodes of 
rapid cyclone development are not common events, suggesting 
that the atmosphere operates in a more subtle manner. 
Specifically, most cases of rapid cyclone development cannot 
be explained using quasi-geostrophy. Therefore, it is 
hypothesized that certain relationships between the afore- 
mentioned physical mechanisms at different spatial and 
temporal scales are crucial for the acceleration of cyclone 
growth and development. Development may also be accelerated 
if both sources of energy (eddy available potential energy 
and zonal kinetic energy) are available to the initial 
perturbation. In other words, a cyclone may develop rapidly 
through a combination of baroclinic and barotropic processes 
if the horizontal and vertical shears are of sufficient 
magnitude. 

This thesis will focus on the processes’ of baroclinic 
and barotropic instabilities by studying the genesis, growth 
and behavior of the Carolinas storm of March 1984. Based 
upon the synoptic analysis and discussion of this cyclone 
(Chapter 3) it is hypothesized that both baroclinic and 
barotropic instabilities contributed to this particular 
cyclone event. Such a hypothesis will be tested in Chapter 


5 using a linear stability analysis. 
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III. SYNOPTIC ANALYSIS AND DISCUSSION 


A. DATA AND METHODS OF ANALYSIS 

Operational height and temperature analyses at 850, 700, 
500, 300 and 200 mb at the 0000 and 1200 GMT synoptic times 
and 3-hourly surface pressure analyses from NMC were used in 
this portion of the study. Height, temperature and wind 
analyses for the same levels from NORAPS were also used. 
Other data include hourly surface observations from the 
first-order stations and ship and buoy reports along coastal 
areas between 282100GMT through 290300GMT. Rawinsonde 
observations from stations throughout the central and 
eastern sections of the USA were also available. 

The ‘hourly observations were manually analyzed for 
temperatures, pressures and pressure tendencies. These anal- 
yses provide a more detailed depiction of the synoptic as 
well as the sub-synoptic events that occurred prior to and 
during genesis, as will be seen in section C of this 
Chapter. 

To obtain information about the vertical structure of 
the atmosphere and the incipient cyclone, four cross- 
sections were chosen. The cross-sections were constructed 
using a technique similar to that described by Whittaker and 
Petersen (1975). Potential temperatures were first calcu- 
lated from the pressure and temperature observations at each 
of the rawinsonde stations along the cross-section. Any 
Static instability at each station iS removed prior to the 
analysis. Temperatures and potential temperatures were then 
linearly interpolated onto 50 vertical levels and 50 hori- 
zontal grid points (all equally spaced). 

The wind cross-section was constructed by PAS E 
computing the wind component perpendicular to the cross- 


Section. These normal components were then linearly 
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interpolated in the vertical for those levels in which wind 
observations were available. Otherwise, the thermal wind 
relation was used to provide an estimate of the winds in the 
following manner. The horizontal temperature gradient was 
first computed using the interpolated temperatures at the 
“two grid points (one on each side) nearest to the station. 
The wind at a level above the highest level in which the 
wind observation was available can then be estimated with 
the assumption that the wind observation at the highest 


available level could be approximated by the geostrophic 


wind. This assumption was verified by geostrophic wind 
calculations (not shown) using the temperature field 
described above. This procedure is repeated for all the 
other levels above. After the winds have been computed to 


100 mb, they are linearly interpolated in the horizontal to 
obtain a 50x50 grid of winds. 

To gain further insight into the structure of the inci- 
pient low as well as possible physical processes, computa- 
tions of static stability, absolute vorticity and Richardson 
numbers were made for each cross-section at 290000GMT. 


These are defined as follows: 


static stability (SS)= -(X/0)(90/3P) (3 


Richardson number = SS/]3U/IP]? (352) 


where oC is the specific volume, 8 the potential tempera- 


ture, P the pressure and U the zonal wind. 


B. SYNOPTIC OVERVIEW 
The 28-30 March 1984 case of rapid cyclone development, 
which produced severe weather and heavy snowfall over the 


eastern portions of the United States, was one of the most 
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destructive on record. A series of tornadoes across South 
and North Carolina caused 57 deaths and 1,248 injuries, 
according to the 1984 Storm Data statistics (Kocin et al., 
1984). Wind-driven tidal surges caused widespread damage 
along the middle Atlantic to the New England coastline. As 
much as 75 cm of snow fell across interior sections of 
Pennsylvania, New York and New England. 

The following discussion will focus on the tracks and 
intensities of the individual low pressure systems that 


existed over a 48-h period between 280000GMT through 


300000GMT. (All future references to time are with respect 
to GMT.) This section describes the synoptic background 
from which the Carolinas storm evolved. A more detailed 


analysis of the genesis period is contained in the next 
section. 

At 280000, a mid- to upper-level synoptic-scaled trough 
extended through the central portions of the US. This 
trough was moving rapidly eastward with little change in 
amplitude until 281200, when it decelerated and began to 
deepen more rapidly (Fig. 3.1). A shortwave trough was also 
observed to the westnorthwest of the major 500 mb low center 
at 281200. The jet streak at 300 mb had winds of 100-140 kt 
and was located over central Texas. 

Three separate low pressure centers can be identified 
from the NMC surface analyses between 280000 and 291200. 
The tracks and intensities of these systems are shown: in 
Figs. 3.2 and 3.3 respectively. At 280000, low I was over 
east Texas with a pressure of 989 mb. This feature was 
imbedded in a large region of low sea-level pressure that 
encompassed most of the midwest sections of the US. Low I 
moved rapidly northeastward while slowly deepening to 984 
mb. At 281200, low II formed over western Tennessee and 
moved slowly eastward to approximately 200 km west of the 
Appalachians. Low II continued to deepen slowly, with the 


central pressure reaching 982 mb 12 h later. 
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Between 281800 through 290000, surface pressure falls 
exceeded 6 mb over a 3-h interval across Georgia and Alabama 
concurrent with the formation and rapid movement of low III. 
Warm, humid air (dew- point depressions less than 3° C) at 
281800 covered the southeastern US as far north as North 
Carolina. Moderate rain was also widespread across the 
eastern half of the US, with snow across New York and 
Pennsylvania. The widespread region of clouds that produced 
the moderate to heavy precipitation can be seen in Fig. 3.4. 

A six-hourly sequence of surface analyses between 281800 
and 290600 shown in Fig. 3.5 vividly depicts the large, 
intensifying cyclonic circulation that occupied the eastern 
half of the United States during this time. Imbedded within 
this large circulation are the three low pressure systems 
(lows I, II and III) described above. The unmarked low was 
a weak, shallow disturbance that had moved eastward across 
North Carolina. . 

Low III developed at 281800 over eastern Alabama along a 
cold front, and moved rapidly northeastward to a position 
over central Georgia by 282100. The available observations 
did not suggest low III originated in the Gulf of Mexico, 
which is consistent with the conclusion in related papers on 
this storm (Kocin et al., 1984; Hillger et al., 1985). A 
very large circular cloud mass southeast of the Appalachians 
can be observed at 282100 from infrared GOES imagery 
(Fig. 3.4). This feature was similar in shape and size’to 
the mesoscale convective complex (MCC) described by Maddox 
(1980). At 282300, the satellite signature of this convec- 
tive area appeared to change shape, perhaps due to a mid or 
upper-level wind maximum. The changes continued until by 
290030 (not shown), a well-defined comma cloud was evident. 
The strong height gradient over this region on the 290000 
NORAPS 500 mb analyses suggests that a wind maximum did 


exist over and to the southwest of the comma cloud 
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Mig. 3.6). Satellite imageries at 282300 and 290100 show 
the area of severe weather near the South Carolina-Georgia 
border (Fig. 3.4). By 290000, the NMC 500 mb analysis (not 
shown) reveals that the shortwave had progressed to a posi- 
tion east of the longwave trough position at 500 mb, and lay 


southwest but in close proximity to an area of severe 


weather. Low I had become primarily a mid-level feature, 
and appeared to be associated with snowfall over 
Pennsylvania, New York and New England at this time. The 


rapid northeastward movement of low III on both Figs. 3.2 
and 3.5 18 apparent until 290000 when its translation speed 
was more consistent with that of a midlatitude synoptic- 
scaled cyclone. 

The size of low III at this time was approximated by 
measuring the diameter of the 984 mb isobar (Fig. 3.5), 


which roughly coincides with the perimeter of the comma- 


shaped cloud (Fig. 3.4). Since the 984 mb isobar does not 
encircle low III, this approximation was only valid to the 
north, east and south of low III. (The extension of this 


isobar to the west of low III reflects the presence of low 
II.) The diameter was measured to be about 1000 km which at 
this latitude approximately corresponds to a disturbance of 
zonal wavenumber 15. 

Between 290300 and 290600, low III deepened very 
rapidly. By 290600, low III became the dominant synoptic 
feature with respect to the other low pressure systems 
discussed above (Fig. 3.5). By noting the change in the 
distance of the 980 mb isobar from the center of low III at 
290000 and 290600 (Fig. 3.5), it was estimated that the 
cyclone doubled in size within this time interval, which 
corresponds to a growth rate of approximately 3/day. The 
980 mb isobar is used in estimating the size change because 
it better depicts the evolution of low III than does the 984 
mb isobar (Big. 3.5)- This growth rate is further 
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substantiated by observing that at 290000, the surface 
pressure perturbation along a 2000 km north-south line 
centered in the middle of the low was approximately -6 mb. 
This pressure perturbation at 290600 became approximately 
-1l mb. These estimates of the growth rates will be 
discussed further in Chapter 5 when compared with the 
results of the linear stability model. 

Between 290300 through 291500, the cyclone (low III) 
deepened from 982 to 965 mb (Fig. 3.3) and grew in size from 
approximately 1000 to 4000 km, while tracking to within 
approximately 100 km northeast of Norfolk, Va. at 290900. 
The most rapid development of the cyclone took place over 
land, although it reached its lowest pressure (at 291500) 
over water. 

Satellite pictures (Fig. 3.4) show that as the cyclone 
developed, the associated convection that was located mostly 
northeastward of low III grew rapidly in horizontal extent. 
The 291200 operational analyses indicate that the cyclone 
extended from the surface through 500 mb. The cyclone was 
now situated completely in the cold air, and had become 
vertically stacked. No phase shift between the 500 mb 
thermal field and height field can be discerned, as had been 
evident at 290000 (Figs. 3.6 and 3.7). Subsequent satellite 
pictures after 291200 (not shown) indicate that the main 
areas of convection were 550.km to the north, east and 
southeast of the cyclone and had a north-south extent: of 
approximately 2000 km. The low began to weaken at 292100, 
while moving northeastward parallel to the coastline. 

The most noteworthy synoptic features at 290000 are the 
distinct westward vertical tilt of the pressure trough , as 
well as the phase difference between the thermal and height 
fields (Figs. 3.5, 3. 6 and omn. These features imply the 
strong baroclinic nature of the atmosphere prior to and 


during genesis. The vertical structure also suggests that 
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the incipient cyclone was primarily a surface disturbance. 
The 850 mb height analysis (Fig. 3.7) further indicates the 
presence of a low-level jet over southeastern Georgia and 
South Carolina and the strong warm advection to the north- 
east of low III. This jet suggests that the low-level 
vertical and horizontal shears were greatly increased. 
Strong vertical shears could result in regions near the 
surface where low Richardson numbers might exist. These 
observations point to shortwave baroclinic instability as a 
possible contributor to the genesis of low III. The strong 
horizontal shears also suggest that barotropic instability 
may have played a role during the early stages of the devel- 
opment of the cyclone. 

The 500 mb wind and temperature fields at 290000 
(Be 3.8) indicate that cold advection is occurring’ at the 
exit region and warm advection at the entrance region of the 
jet. This suggests that rising motion would occur under- 
neath and sinking motion to the south of-the exit region of 
the 500 mb jet (Shapiro, 1982). In fact, most of. the severe 
convection occurred underneath and to the north side of this 
region. Low III was also located almost directly under this 
area. 

In summary, the Carolinas storm appeared to have formed 
in eastern Alabama at 281800 (low III in Fig. 3.2). The 
cyclone moved rapidly northeastward until 290000, when it 
began to decelerate quickly. The low continued to move‘at 
this slower speed until its demise. Furthermore, the 
cyclone began to deepen rapidly after 290000. Most of the 
Severe convection occurred between 290000 and 290600. 
However, the convection was confined to a small region 
mostly to the northeast: of the cyclone. The atmosphere 
appeared to be highly baroclinic at 290000 so that baro- 
clinic processes could be important in the formation of low 


III and perhaps its subsequent rapid development. This idea 
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will be explored further using a linear stability model 
(Chapter 5). It is also apparent that to understand the 
behavior of this particular case of explosive cyclone devel- 
opment, a more detailed analysis of the environment near 
290000 is warranted. 


C.. DETAILED SURFACE ANALYSES OF THE GENESIS PERIOD 

To provide a more coherent synopsis of the events that 
occurred around 290000, hourly surface observations in the 
eastern US for seven time periods (282100 to 290300) were 
analyzed. Each analysis includes the temperatures, pres- 
sures, pressure tendencies and frontal positions. These 
analyses augment Fig. 3.5 and were crucially aided by the 
information gleaned from satellite pictures. 

The 282100 analysis (Fig. 3.9) shows a strong tempera- 
ture gradient oriented along the lee slopes of the 
Appalachians. Low III was imbedded in this gradient over 
northeastern Alabama, with temperatures greater than 80° F 
to the southeast and less than 60° F to the northwest. The 
temperature analysis also suggests an extensive amount of 
cold air damming, which might have enhanced the temperature 
gradient. At this time, very little warm advection appears 
to be taking place ahead of the low. 

Surface damming of cold air east of the Appalachians is 
particularly apparent north of the warm front. Richwien 
(1980) found that damming frequently occurs after a cold 
anticyclone becomes established over the northeastern US. 
The ridge axis typically extends southward between the 
Appalachians and the east coast. In this case, the surface 
analysis shows a weaker ridge across Virginia than Richwien 
observed, extending from a moderately strong anticyclone 
situated well to the north of the Great Lakes. This weak 
ridge over Virginia is coincident with colder air to the 
north of low III (Fig. 3.9). Consequently, it seems reason- 


able that damming may also have contributed to the rapid 
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development of low III. However, damming may also have 
played an important role in the formation, location and 
amount of snow to the north. 


Detmeens 252000 and 282100, rapid pressure falls (rises) 


ahead of (behind) low III can be readily observed 
(Fig. 3.10), which was largely related to the movement of 
this cyclone. In fact, the movement of low III is almost 


parallel to the isallobaric vector, which is directed north- 
eastward from north-central Georgia to. northema South 
Carolina. 

Major changes in the surface pressure and temperature 
fields appeared to. have taken place between 282100 and 
232200 (Fig. 3.11). Low III had moved very rapidly east- 
northeast to a position over central South Carolina. The 
central pressure fell to 978 mb and the strong baroclinic 
zone to the north of the low was perturbed. Also of partic- 
ularly interest in Fig. 3.11 is the warm air now prevalent 
ahead of the low. Large pressure falls are found to the 
aig@eheast and™south of low ITI. (Fig. 3.12). The isallobaric 
vector is now directed east-northeast from northeastern 
Georgia towards the North/South Carolina border, which again 
is almost parallel to the movement of low III. The movement 
of low III appears to be towards the region of largest 
hourly surface pressure falls (Fig. 3.12). 

By 282300, low III was located over northeastern South 


Carolina, continuing its rapid movement over the past hour 
ie. 3.13). Furthermore, the central pressure had 
decreased an additional 2 mb. Warm air advection had become 


much more prevalent ahead of the low, but the baroclinic 
zone was less perturbed than in the previous hour. One 
possible explanation for this is the apparent cold dome 
associated with the downdraft at Greensboro, North Carolina 
where moderate rain was reported at this time. The tempera- 


ture has decreased by 2° F over the past hour at this 
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station. This drop in temperature had helped to strengthen 
the baroclinic zone to the north of the low. Pressure falls 
are again noted over a large area to the north, east and 
south of the low (Fig. 3.14) in conjunction with the move- 
ment and development of this feature. However, the isallo- 
baric vector is not well-defined. The movement of low III 
is also not directly towards well-defined regions of pres- 
sure falls. 

The 290000 surface analysis (Fig. 3.5) indicates that 
the low was situated at the North/South Carolina border, 
giving a forward speed of about 30 kt. Subsequent analyses 
show that low III moved at approximately this speed. until 
its demise. Warm advection was still prevalent at 290000 
ahead of the low, and colder air had begun to move southward 
toward the low from the north (not shown). The baroclinic 
zone also gradually weakened with time. The pressure anal- 
ysis at 290000 (Fig. 3.5) suggests that the low filled 
slightly, but the surface pressure tendencies between 0000 
and 0100 indicate large pressure falls northeast of the low 
(not shown), corresponding to the northeastward movement of 
the low. . 
| These surface analyses in conjunction with Fig. 3.4 have 
shown that between 281800-282200, low III was primarily 
located underneath its associated cloud mass. At 282300, 
the low became removed from this cloud mass. Furthermore, 
the satellite signature at 282300 became more typical of a 
frontal wave, suggesting that baroclinic instability might 
be important. It appears that the baroclinity of the atmos- 
phere in this region at this time is increased due to the 
presence of both the upper-level jet streak and the low- 
level jet. The horizontal and vertical shears associated 
with these jets together with the convectively unstable 
atmosphere appear to have created a favorable environment 
for the low to extract energy from this enhanced 


baroclinity. 
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D. VERTICAL CROSS-SECTIONS 

To understand better the stability characteristics, 
structure of frontal zones and wind and temperature profiles 
associated with the 1984 Carolinas storm, four cross- 
sections (Fig. 3.15) were constructed prior to and following 
the genesis of low III. The first and second cross-sections 
were chosen perpendicular to the positions of the 500 mb jet 
at 281200 and 290000 respectively. The third cross-section 
is closely aligned with the upper-level trough, while the 
fourth goes through low III and the areas of severe convec- 
tion at 290000. 

The first cross-section is constructed at 281200 between 
Salem, Illinois and Brownsville, Texas to provide insights 
into the state of the atmosphere prior to genesis. However, 
analyses from this cross-section do not contribute signifi- 
cantly to the discussion at 290000, and therefore will not 
be presented. However it is noted that the upper-level jet 
streak was located over central Texas exhibiting a similar 
intensity to that observed at 290000. It is therefore 
concluded that the jet streak did not undergo any signifi- 
cant changes in intensity between 281200 and 290000. 

The second cross-section (Fig. ss. 169 from Peoria, 
Illinois (PIA) and West Palm Beach, Florida (PBI) shows an 
area of very strong vertical shear near the surface between 
PBI and Athens, Georgia (AHN). There is also the suggestion 
of a low-level jet at 750 mb between AYS and PBI. A well- 
defined jet is found near AHN at approximately 250 mb. The 
surface cold front can be seen north of Waycross, Georgia 
CABS). This position is consistent with its analyzed posi- 
tion on Fig. 3.5 at 290000. 

Low static stability values shown in Fig. 3.17 south of 
AHN correspond well with the areas of convection at this 
time. Note that areas of low static stability near the 


surface existed north and south of the cold front. 
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Richardson numbers (Ri) depicted in Fig. 3.18 show regions 
of low Ri at 800 mb over Athens and to the south near the 
surface in association with the low-level jet. Other areas 
of low Ri exist beneath the upper-level jet. These low 
values of Ri would therefore suggest that the environment 
was favorable for the genesis of a subsynoptic or meso-alpha 
scaled disturbance. The regions of lower Ri values near the 
surface should also favor baroclinic growth. In fact. To 
III was located approximately 250 km to the northeast of 
this cross-section between AHN and AYS at this time but was 
almost directly over Athens 2 hours previous. Absolute 
vorticity values over a large area south of AYS above 700 mb 
are small (Fie. 3.1992 This area is on the anticyclonic 
side of the very strong jet over AHN. Inertial instability 
might therefore have been present. 

The third cross-section (Fig. 3.20) was chosen to 
provide a better definition of the location of the jet 
streak over the southeast US by comparing it with profile II 
(Fig: 3.16). This cross-section from Dayton, Ohio (DAY) to 
Brownsville, Texas (BRO) shows a jet at Jackson, Mississippi 
(JAN) at around 250 mb. Note that in this profile the jet 
is weaker than the one found in profile II (Fig. 3.16) and 
located at a higher elevation. A lower-level jet (near 400 
mb) appears over Lake Charles, Louisiana (LCH) at 290000. 
The strongest surface temperature gradient exists between 
LCH and BRO. 

Low static stability values near the surface south of 


BNA (not shown) are indicative of the push of cold air 


southward behind the cold front. Near the surface south of 
JAN, values of Ri are relatively small (not shown). Areas 
of low Ri values are found at 400 mb near JAN. This is a 


result of vertical shear underneath and to the north of the 
jet. Absolute vorticity values south of LCH in Fig. 3.21 


show that the jet at this location is more inertially stable 
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Cam that over AYS (Fig. 3.19). Therefore, the jet closer 
to the area where low III formed is more inertially 
unstable, providing further evidence that the jet streak was 
located in profile II. 


The fourth cross-section (Eig: 3.22) extends from 
Appalachicola, Florida (AQQ) to Chatham, Massachusetts 
(CHH). In this profile, a jet is well-defined over Cape 
Hatteras (HAT) at about 200 mb, as well as a lower-level jet 
around 500 mb. A strong warm front can be seen between HAT 
and Wallops, Virginia (WAL). Low III existed between Cape 
Hatteras and Charleston, South Carolina (CHS), as can be 


seen from the isotach pattern. 


Along this profile, regions of low static stability 


(Fig. 3.23) are evident near the surface. Severe convection 
was occurring around HAT at this time. Low Ri (not shown) 
prevailed over WAL. Absolute vorticity values were very 


small between HAT and CHS; particularly at 250 mb 
(Fig. 3.24), which corresponds with that area closest to the 
location of low III. 

In summary, the atmosphere at 290000 was characterized 
by strong jets and associated baroclinic zones, areas of low 
Ri and two regions where inertial instability may have 
existed. In addition, the horizontal extent of the cyclone 
was approximately equivalent to that of a wave disturbance 
of zonal wavenumber 15 and was mostly confined to the low 
levels. Gall (1976a) has shown that for favorable linear 
growth of sub-synoptic scale baroclinic wave disturbances, 
the zonal Ri should be smaller near the surface than aloft. 
Mullen (1979) also found from 22 cases of polar lows that 
the mean center of the surface cyclone during genesis was 
characterized by Ri on the order of 1. Therefore, shortwave 
baroclinic instability may have played an important role in 
the genesis of this cyclone. The presence of the strong jet 


streak greatly increased the horizontal and vertical wind 
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shears near the location of the cyclone. As a result, baro- 
tropic instability may also have contributed to the genesis 
of low III. These ideas will be further tested in Chapter 5. 


E. NUMERICAL MODEL PERFORMANCE 

Previous studies have shown that numerical forecasts 
usually are not very successful in cases of rapid cyclone 
development (Sanders and Gyakum,1980; and others). Lites 
therefore of interest to evaluate the performance of the 
operational numerical models for this particular cyclone 
event. 

The 36-h surface pressure forecasts from the NOGAPS 
(Rosmond, 1981), NORAPS (Hodur, 1982) and the NMC LFM were 
compared with those from the NMC surface analyses 
uid Cae oul All of the models predicted a general decrease 
in surface pressure of the main low pressure system, but 
failed to deepen the low to the extent observed. The fore- 
casts of the speed of movement were also slow (not shown). 
In addition, not all of these models predicted the genesis 
of low III. Only NORAPS and MASS (Kocin et al., 1984) 
correctly predicted such an event. Kocin et al. (1984) 
suggested that incorrect specification of the vertical 
latent heat distribution and wind fields may have been 
responsible for some of the forecast errors in these models. 
NOGAPS forecasts slightly overdeveloped low III when it had 
reached the mature stage over water, but was too weak when 
the incipient low was still over land (not shown). These 
results are consistent with those obtained by Toll and Clune 
(1985) who found that NOGAPS forecasts of the surface pres- 
sure of maritime low centers are too weak during the early 
stages of development, but are too strong during the mature 
and dissipating stages. 

A more thorough examination of NOGAPS was conducted 
using a vortex tracking program developed by Williamson 
(1981) and modified by Brody et al. (1984) and Harr and 
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sua (1985). The cyclone positions forecast by NOGAPS are 
primarily to the left and behind the analyzed positions 
(Fig. 3.25). The slower-than-observed movement of the model 
cyclone is better depicted in Fig. 3.26. These results 
imply that NOGAPS was slow in predicting the foreward move- 
ment of the storm. 

These observations imply that the regional numerical 
models appear to be better. capable of predicting the 


behavior of a cyclone prior to and during genesis than the 


global models. This might be anticipated in view of the 
fact that the incipient cyclone in this case was on the 
meso-alpha scale. In fact, NORAPS (Hodur, personal communi- 


cation) produced more accurate forecasts of both the posi- 
tion and intensity of low III when the model was initialized 
with a high resolution analysis and its 12-h forecasts from 
the previous 12 h. (The NORAPS forecasts evaluated in this 
study were initialized by NOGAPS.) The more accurate NORAPS 
forecasts will be used in Chapter 5 as part of a linear 


stability study of this cyclone. 
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Ele cs 281200 NORAPS o mb analysis of heights 
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Fig. 





3.4 Selected SMS-GOES infrared satellite imagery 


between 281800 and 290330. The surface low posi- 


tions are indicated by L and c refers to the 
expanding region of lower cloud top temperatures. 
(Reproduced From Kocin et al., 1984). 
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Pipes. 5 Sea-level pressure (mb) and surface frontal 
analyses at 281800, 290000 and 290600. Three low 
es centers indicated on as 3.2 are also 

abelled. (Reproduced from Kocin eť al., 1984). 
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Pig. 3.8 Isotach (solid, kt) and temperature (dashed 
o C) analyses at 500 mb at 290000. The letter L 
refers to the center of the 500 mb low center. 
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PLENAS Surface analysis for the southeastern United 
States at 282100. Solid lines_are isobars (mb) and 
dashed lines are isotherms (° F). 
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Fig. 





aS Horizontal orientations and extent of the 
four cross-sections used pia this study. 
Cross-section I includes rawinsonde data from. 
Brownsville (Texas), Victoria (Texas). Longview 
Texas), Little. Rock (Arkansas ) and. alem 
lliinois) as for 2c 200F Cross-section II_includes 
rawinsonde data from West Palm Beach a 
aycross (Georgia), Athens (Georgia), Nashville 
Tennessee), Salem (Illinois) and Peorla (Illinois) 
or 290000. Cross-section [II includes rawinsonde 
data from Brownsville Texas), Lake Charles 
Louisiana), Jackson (Mississippi) Nashville 
Tennessee and Dayton (Green for 290000. 
ross-section IV includes rawinsonde data from 
Appalachicola (Florida, Waycross (Georgia), 
Charleston (South Carolina), Cape Hatteras (Nor 
Carolina), Wallops Island (Virginia) and Chatham 
(Massachusetts) for 290000. 
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Fig. 3.17 Static Stability (197º m? kg! Pa! ) for 
the cross-section in Fig. 3.16. 
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Fig. 3.19 Absolute vorticity (10% s"" ) for the 
cross-section in Fig. 3.16. 
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IV. HYPOTHESIS 


In this case of rapid cyclone development over the 
Carolinas, a strongly baroclinic environment was observed 
prior to genesis. Furthermore, an area of low static 
stability was present downstream of the polar front jet. 
The existence of the jet streak increased the horizontal and 
vertical shears. Therefore, both baroclinic and barotropic 
instabilities may have been important in the genesis of this 
cyclone. As the cyclone went through its life cycle, 
convective processes appeared to have played an increasingly 
important role in the development and maintenance process of 
the cyclone. Although the Carolinas storm was over land, it 
exhibited many characteristics of a rapidly deepening mari- 
time cyclone. Therefore, it 1s reasonable that these mecha- 
nisms may also play similar roles for cyclones over water. 
Based on the analyses in Chapter 3, the following sequence 
of events iS proposed to explain the genesis and rapid 
development of maritime cyclones. 

An intense upper-level jet in thermal balance with a 
strong low-level baroclinic environment first moves into a 
convectively unstable region. Palmen and Newton (1969) 
showed that a strong divergence field is likely to be found 
at the left exit region of the jet. The movement of the jet 
also represents a rapid increase of kinetic energy in the 
region. As the jet streak moves sto a posiciones 
Synoptic scale trough, the strong divergence associated with 


the left front quadrant of the jet is augmented by the 


ambient divergence of the trough. The enhanced divergence 
pattern in this quadrant induces a surface cyclone. 
Concurrently, the import of kinetic energy invigorates the 


atmosphere in such a way as to initiate and maintain the 
convection, particularly if the atmosphere is convectively 
unstable. In other words, the atmosphere attempts to reduce 


the strong vertical wind shear associated with the jet via 
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cenveetion. Strong vertical wind shear and low static 
stability reduce the Richardson number so that shortwave 
baroclinic instability can become an important contributor 
to the development process. The characteristic spatial 
scales of this mechanism suggest that the surface cyclone 
will form on. the meso-alpha, rather than the synoptic, 
scale. On the meso-alpha scale, a strongly baroclinic envi- 


ronment is able to organize the convective processes more 


effectively. Furthermore, the baroclinic energy. conversions 
can be better realized. As a result, convection is effi- 
ciently organized around the incipient cyclone. The strong 


Shears associated with the jet streak can also create a 
barotropically unstable atmosphere in both the vertical and 
horizontal. Therefore, the growing perturbation can extract 
eddy available potential energy through baroclinic processes 
and zonal kinetic energy through barotropic mechanisms. The 
large amount of latent heat release near the incipient 
cyclone provides strong diabatic heating, which can 
partially offset the self-limiting process during the devel- 
opment phase (Palmen and Newton, 1969). 

As the perturbation reaches its mature stage, the 
cyclone becomes more vertically aligned, which Suggests that 
baroclinic instability has become a less dominant factor. 
Rapid development slows when convection decreases in inten- 
sity, or when the main areas of convection move far enough 
away from the cyclone as to weaken the effect of diabatic 
heating. The maturing cyclone is typically in the cold air 
with a mid-level trough situated above. Thus, according to 
the Petterssen development equation (Petterssen, 1956), the 
cyclone would quickly cease its development and begin to 
weaken. During convection, both the vertical and horizontal 
Shears are weakened, thus reducing the possibility of baro- 


Eap 1 Geanstabl lity, 
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In conclusion, genesis and rapid development is initi- 
ated by strong baroclinic instability. This instability 
favors the existence of a strong jet streak. However, a 
very intense jet streak can also produce barotropic insta- 
bility through a transfer of momentum from the mean flow to 
the eddy in both the horizontal and vertical. Baroclinic 
and barotropic instabilities may therefore combine in such a 


way as to accelerate the growth process by providing the 


incipient cyclone with two energy sources. This hybrid 
perturbation can also be further assisted by moist 
processes. 


In the particular case of the Carolinas storm, NORAPS 
analyses at 290000 show the incipient cyclone to be below 
850 mb (Figs. 3460 Anes At 500 mb, a jet streak (110 
kt) had moved around the base of the trough at this time, 
and was approximately 350 km southwest of the surface 
cyclone position. The 850 mb and 700 mb temperature and 
moisture analyses (not shown) as well as static stability 
calculations (Fig.3.17) reveal that the atmosphere over most 
of the area around the surface cyclone was’ convectively 
unstable. A low-level jet and strong warm air advection at 
850 mb were found to extend to the east of the surface 
cyclone. The distinctive lag between the 500 mb isotherm 
and height pattern of approximately 45 degrees at 281200 had 
decreased to about 20-30 degrees by 290000. Furthermore, 
the 500 mb trough had significantly deepened over this same 
time interval. The cross-section analyses indicate that 
strong vertical shears below 850 mb existed near the area of 
genesis and rapid development, where low values of the 
Richardson number prevailed (Fig.3.18). Low values of abso- 
lute vorticity (Fig.3.19) were observed at mid and upper 
levels on the anticyclonic side of the jet. 

Therefore, all of these conditions are consistent with 


the above proposal so that baroclinic and barotropic 
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instabilities appear to be responsible for the genesis of 
the Carolinas storm at 290000. The main purpose of this 
thesis 1S to examine such a hypothesis and to determine the 
role and relative importance of shortwave baroclinic and 
barotropic instabilities in the genesis and rapid 
development of this cyclone. Gall (1976a) has shown how 
these two instabilities can be diagnosed with a linear 
stability model. À version of this model will therefore be 
used to determine whether these mechanisms can help explain 


the behavior of this cyclone. 
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V. LINEAR STABILITY ANALYSIS 


A. INTRODUCTION | 

The analyses in Chapter 3 led to the hypothesis in 
Chapter 4 that the processes of baroclinic and barotropic 
instabilities might explain the rapid growth and development 
of the Carolinas storm. Specifically, development was 
accelerated by two energy conversion processes operating in 
tandem from two different Sources. These ideas will be 
tested using a linear stability model. Different versions 
of the model have been described and used by Gall 
(1976a,b,c), Blakeslee and Gall 1978) and Gall et al. 
(19379). 

Linear theory considers the perturbations to be infini- 
tesimally small, while observed eddies in the atmosphere are 
typically of large amplitude. Thus, lt is’ assumed that 
disturbances, as they grow from small to observable size, do 
not appreciably change their structure. Hence, linear 
stability theory can be extended to those perturbations of 
finite amplitude. In fact, Gall (1976a) suggested that a 
linear model can produce perturbations with structures 
Similar to those observed inthe atmosphere. The linear 
stability analysis should provide a reasonable estimate of 
the growth rates, phase speeds, structure and energetics of 
a perturbation introduced to an unstable mean state that is 
representative of the atmosphere during the genesis of the 
Carolinas storm. As discussed by Gall (1976a), this tech- 
nique also includes the effect of ageostrophic motions. 
Therefore, linear theory should be applicable to the 
Carolinas storm during the genesis and early stages of 
development. However, this linear analysis cannot provide 
an estimate of the location of the perturbation that would 


grow from the unstable initial state. Such an estimate can 
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be derived from vorticity advection arguments. On the other 
hand, the behavior and characteristics of the perturbation 
can be identified from the model. Therefore, a combination 
of these two approaches might provide better estimates of 
both the position and intensity of the incipient cyclone. 

For the results from the linear analysis to be essen- 
tially valid, it is important , to exclude any significant 
energy source that might be available from a nonlinear 
mechanism. For example, the initial conditions must be 
constructed along a cross-section that minimizes the nonli- 
near effects of moisture and convective processes. Itris 
also possible that nonlinearities or unstable modes other 
than the desired Rossby modes would dominate the growth 
Spectrum, and thus be of little use for the purposes of this 
study. Therefore, extra precautions become necessary to 
ensure that only the unstable baroclinic and barotropic 
modes emerge after the requisite number of integrations have 
been completed. 

In the next section, the linear stability model.to be 
used will be described. In section C, the cross-sections of 
zonal wind and temperature to be studied are presented. The 
results of the linear analysis, including the growth rates, 
energetics and structure of the perturbations, are shown in 


section D and compared with the observed fields. 


B. DESCRIPTION OF THE LINEAR STABILITY MODEL . 
The linear stability model provided by Professor R. L. 
Gl was first developed by him in 1976 (Gall, 1976a). 
The model is structured around the linearized primitive 
equations in pressure coordinates, which retain the effects 
of spherical geometry. The equations are written in 'semi- 


spectral" form (Saltzman, 1957), in which the perturbations 
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to the mean zonal flow are expressed as a Fourier series in 
zonal wavenumber. The meridional and vertical derivatives 
are evaluated by an energy-conserving centered difference 
scheme, while the zonal derivatives are evaluated analyt- 
ically (consistent with the assumption that the atmosphere 
is invariant in the zonal direction). The time derivatives 
are replaced by centered differences with time-smoothing 
(Bleck, 1974) applied at each time step to reduce high 
frequency oscillations due to gravity waves, and to suppress 
the computational mode in time. The velocity components are 
staggered on the grid to reduce the spatial computational 
mode. 

The model is confined to a channel in the meridional 
direction. No mass flux through the channel walls is 
allowed, and the frictional stresses and the turbulent 
transfer of heat and moisture across the walls are zero. No 
topography is included on the lower boundary. At the 
surface, only the vertical turbulent transfer of momentum is 
included, but none for heat or moisture. Thus, the total 
energy within the channel slowly diminishes during the 
course of the experiment. 

The two dimensional (y,p) grid is oriented perpendicular 
to the zonal wind. To initiate the development of wave 
motion, a temperature perturbation of -0.1 to +0.1° Kelvin 
1s added to the zonal mean temperature. After each 
timestep, the flow is modified by restoring the zonal flow 
to its original values and restricting the eddy portion of 
the flow to one zonal wavenumber. Therefore, the relative 
growth rates of different waves passing through the grid can 
be calculated. The growth rate of each wave then gives an 
estimate of the amount of the instability of the zonal flow. 

The method for obtaining the linear solutions is the 
initial value technique. Specifically, an initial guess for 


the Fourier coefficients at each grid point of the 
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north-south cross-section is made, and the equations for 
these coefficients are numerically integrated in time until 
the most unstable mode dominates the flow, and the real and 


imaginary parts of the phase speed are the same at all grid 


points. 
With the linear primitive equations, gravity-wave type 
solutions must be addressed. The introduction of a wave in 


the temperature field represents an imbalance to the geos- 
trophic balance imposed on the mean zonal flow. As a 
result, gravity waves will be generated early in the experi- 
ment. These waves will remain in the solutions throughout 
the integration because of the absence of frictional forces 
or interactions with the zonal flow. The width of the jet 
determines the scale of the perturbation. Consequently, in 
other regions of the grid, the amplitude of the gravity 
waves may exceed that of the unstable wave. Thus, equality 
of the phase speeds and growth rates at all grid points 
cannot be expected. For these reasons, the integration must 
continue until the amplitude of the Rossby wave sufficiently 
exceeds that of the gravity waves and all other more slowly 
growing Rossby modes over most of the grid. Ideally, both 
growth rate and phase speed should become constant after 
development of the wave becomes exponential. 

It should be noted that if two waves with the same zonal 
wavenumber have very similar growth rates, the initial value 
technique will be unable to discriminate between them. -In 
this case, the mode with the most rapid growth rate will 
always be the larger of the two and will dominate the flow 
if given enough time. Thus, the model can only isolate the 
most unstable mode. In other words, other unstable modes 
cannot be identified by this technique even though they may 
exist. However, the emergence of this most unstable mode 
may not occur until the perturbation amplitude reaches 


unrealistically large values and the numerical assumptions 


6/ 


break down. For instance, the surface pressure perturba- 
tions could become large enough to give negative total 
pressure. 

In this study, Call sao model was modified such 
that integrations using various horizontal and vertical 
spatial resolutions can be made with minimal difficulty. 
Depending on the initial conditions, the model grid for 
these experiments consists of 10 or ll grid points in both 
the horizontal and the vertical. This resolution was neces- 
sary to resolve a channel that extends from 26-40° N 
(approximately 1800 km). Thus, the horizontal grid length 
is approximately 180 km, compared with the 200 km used in 
Gall (1976c). Various horizontal and vertical grid lengths 
have been tested and no significant differences can be iden- 
tified. Therefore, results using the above grid resolution 
for all experiments will be presented. The relatively 
narrow width of the channel used in this study compared with 
that of Gall (1976c) justified the f-plane approximation, 
with the value for the Coriolis force being valid at 33° N. 
Gall's calculations were based on a beta plane, with the 
channel extending from approximately 25-60° N. The more 
southward location of the present model channel for these 
experiments makes the flow more susceptible to inertial 


instability than that in Gall's experiments. 


C. INITIALIZATION 

Four different sets of initial conditions were used in 
this stability analysis. They include the real data, NORAPS 
analysis, NORAPS 24-h forecast and a data set generated from 
analytic functions. The first three sets of initial condi- 
tions were chosen so that the cross-section was oriented 
perpendicular to the jet streak, with the wind maximum 
located in the center of the profile. The cross-sections 
for both NORAPS cases closely approximated the location of 


the cross-section using the real data. It is also important 
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to emphasize that to be consistent with the design of this 
stability analysis model, the cross-sections do not include 
significant changes in terrain or major convective areas. 
The analytical conditions were chosen for sensitivity tests. 

The real data case uses the information from the 50x50 
grid of profile II described in Chapter 3 (Fig.3.16). The 
geostrophic winds were first calculated at each rawinsonde 
station from the horizontal temperature gradient across the 
station using the thermal wind equation. The winds were 
then linearly interpolated between stations. These geos- 
trophic winds were then interpolated horizontally using a 
cubic spline onto the three staggered grids used by the 
model, with a maximum of 10 horizontal grid points and 10 
vertical levels. The resulting horizontal grid length was 
approximately 178 kn, with the horizontal width of the 
cross-section equal to approximately 1780 km. The vertical 
levels were placed at 1000, 909, 787, 715, 590, 511, 402, 
301, 205 and 100 mb. These levels were chosen to eliminate 
- interpolation in the vertical and to approximate the stan- 
dard atmospheric levels. The temperature at each grid point 
was then calculated from these geostrophic winds. In this 
way, the winds and temperatures were kept in geostrophic 
balance. 

To help isolate the Rossby mode, static and inertial 
stabilities were checked and adjusted as necessary such that 
the model atmosphere was both statically and inertiaily 
stable. A one-sided differencing scheme was used to calcu- 
late both the static stability and absolute vorticity at 
each grid point. If the static instability is found at a 
grid point, the potential temperature at the grid point one 
level below was changed such that stability prevailed with 
respect to the potential temperature at the grid point imme- 
diately above. When the absolute vorticity at a grid point 


is negative (that is, inertially unstable), the geostrophic 
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wind was changed at the grid point immediately to the south 
Such that inertial stability existed with respect to the 
geostrophic wind at the original grid point. The wind and 
temperature fields were rebalanced geostrophically. Another 
stability check was then invoked, but no further instabili- 
ties were found. 

The NORAPS analysis at 290000 provided another set of 
initial conditions for profile II in Fig. 3.15. To minimize 
the amount of interpolation necessary, the primary model 
grid points in both the horizontal and vertical were chosen 
to coincide with the NORAPS analyses. In this way, only one 
interpolation in the horizontal and one in the vertical was 
necessary. This results in a horizontal grid length of 179 
km and 11 vertical levels (1000, 925, 850, 700, 500, 400, 
300, 250, 200, 150, and 100 mb). To obtain the initial 
conditions, geostrophic winds were calculated from the 
NORAPS-analyzed temperatures. As in the real data case, a 
check for static and inertial instabilities was invoked for 
this cross-section. However, only changes to establish 
inertial stability were necessary as the NORAPS temperature 
analysis was statically stable. The resultant cross-section 
is shown in Fig. 5.1. l 

A third cross-section was produced using a NORAPS 24-h 
forecast initialized at 280000. This forecast is different 
from the one evaluated in Chapter 3. For the forecast used 
here, the first-guess fields consist of the 12-h NORAPS 
forecast from the previous 12 h. On the other hand, the 
operational NORAPS forecasts discussed in Chapter 3 used a 
6-h NOGAPS forecast as the first-guess field. The former 
analysis also included more observations (Hodur, personal 
communication). The forecasts used in this present analysis 
provided a better prediction of the intensity of low III, 
but the location of the low center was to the southwest of 


the analyzed position. Nevertheless, this predicted 
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location of the cyclone was in better agreement with the 
analysis than that from operational forecasts. The grid 
points in this cross-section are identical to those for the 
NORAPS analysis case. Geostrophic winds were calculated 
from this temperature field (Fig. 5.2). Static and inertial 
stability checks were also invoked. | 

The temperature and wind fields in Figs. 5.1 and 5.2 are 
very Similar, except that the maximum geostrophic winds for 
the NORAPS analysis case are 10-15 m/s stronger than those 
for the forecast. Stronger horizontal shear exists at the 
jet level for the analysis case, whereas the forecast shows 
stronger horizontal shears near the 700-850 mb levels. 
Strong vertical shear can also be seen near AHN (Athens, 
Georgia) in both cases. Note further that in both cases, 
the maximum winds are weaker than the observed zonal winds 
Pi p ows). The positions of the jet in the analysis and 
the forecast are to the north of that observed. However, it 
must be pointed out that the winds in Figs. 5.1 and 5.2 were 
geostrophically computed. Timfact; the geostrophic winds 
derived from the temperature field in Fig. 3.16 have similar 
magnitudes as those in the NORAPS analysis. The differences 
between the forecast and the analysis can be studied further 
by examining the surface charts (Figs 5.3 and 5.4). Note 
that the cross-section goes through the center of the inci- 
pient low in the forecast (Fig. 5.3), but is positioned to 
the southwest of the analyzed cyclone (Fig. 5.4). Therefore, 
the signature of the low-level perturbation is better 
defined in the forecast cross-section, as can be seen from 
the isotachs to the left (northwest) of AHN, where a 
cyclonic horizontal shear exists (Fig. 5.2). 

The fourth set of initial conditions were analytically 
computed in a way similar to that in Gall (1976c). These 
conditions (Fig. 5.5) were used to test the sensitivity of 


the model for different jet intensities. The zonal mean 
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pressure field is independent of latitude, with the surface 
pressure set at 1013 mb. The horizontal grid is identical 
to that of the NORAPS cases. The ten vertical levels are 
from 1000 to 100 mb in increments of 100 mb. The winds are 


computed using the following equations: 


U(y,p)=1/2(U[(p)cos(Trx/L)) (om) 
Uelp)=((1-p/po)/.775)UMAX (DIAZ OA 5 
U e (p )=UMAX (p < 225 mb) (5:31) 


where UMAX is the maximum value of the zonal mean wind 


UCy> pI L the width of the channels ande p themmis@miace 
pressure. In this study, three different values of UMAX 
(40, 60 and 80 m/s) were used. The latter wind maximum was 
chosen to correspond with that of the other initial 
conditions. The temperature field is computed from the 
zonal wind using the thermal wind equation. The boundary 


condition required to solve (5.3) is given by the zonal mean 


temperature T(p) at the center of the channel: 


T(p)=273.2(p/pa) (p > 225 mb) (5.4) 


T(P)=205 378% (p < 225 mb) (Due DD) 


The resultant vertical temperature profile at the center 
of the channel is a linear function of height, with a 
constant lapse rate to approximately 225 mb. Above 225 mb, 


the temperatures increase with height. The zonal wind is a 


72 


maximum close to the center of the channel above 225 mb. 
The zonal mean temperature and wind fields are in geos- 
trophic balance. 

In addition to the static and inertial stability checks 
already discussed, calculations of the absolute MOL TEL y 
gradient were made to diagnose whether the initial condi- 
tions in each of the four data sets were barotropically 
unstable. The 80 m/s jet in the analytic case was iner- 
tially unstable, but barotropically stable. For the other 
three initial conditions, barotropic instability is preva- 
Tent to the south of the jet. The major difference is in 
the elevation at which barotropic instability existed. For 
the real data and NORAPS analysis case, this instability 
existed at the level of the 250 mb jet. However, in the 
NORAPS forecast case this instability was found mostly at 
the 850 mb level south of Athens. As will be discussed in 
felon seer rom, “these ditferences contributed"to Ssigniti-- 
cant differences in the growth rate spectrum. 

The model was integrated for each of the four initial 
conditions for a spectrum of wavenumbers to determine the 
most unstable wave. The structure and energetics of the 
most unstable mode will then be examined. The results from 


these experiments will be discussed in the next section. 


D. RESULTS 
1. Growth Rates and Phase Speed 

The growth rates of the perturbations from the 
analytic case, NORAPS analysis and the NORAPS forecast are 
illustrated in Fig. 5.6. In the real data case, it appears 
that smaller scale instabilities contaminated the spectrum 
so that the growth rates of the most unstable Rossby modes 
became rather unrealistic. The results from the real data 
experiments will therefore not be shown, but will be quali- 


tatively discussed. 
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The growth rate spectra for both NORAPS cases indi- 
cate that the wavelength of the most unstable perturbation 
occurs around wavenumbers 10-20, which is equivalent toa 
wavelength of approximately 2000 km at NS For these 
wavenumbers, the growth rate is about 3.2-3.5 per day, which 
corresponds to a doubling time of approximately .23 days 
(Sion). Both of these growth rate curves exhibit similar 
tendencies, with lowest growth rates occurring at wavenumber 
de In the NORAPS analysis case, growth rates slowly 
decrease to wavenumber 30, compared to more constant growth 
rates in the NORAPS forecast case at higher wavenumbers. 
The growth rates for the analytic case with UMAX equal to 80 
m/s are approximately constant between wavenumbers 3 and 15, 
but start to decrease at higher wavenumbers. The maximum 
growth rate (approximately 2/day) is also much less than 
that for the NORAPS cases. The curve for the analytic case 
changes dramatically for different jet intensities. For a 
40 m/s jet, the growth rate spectrum is very similar to that 
of Gall (1976c) and Hodur (1984), with a sharp peak at 
wavenumber 15 and a secondary maxima at wavenumber 40. a 
the 60 m/s case, the growth rate curve is flatter, although 
a maximum is still at wavenumber 15. As will be seen later, 
these changes in the growth rate spectrum may be related to 
the energetics. In the weaker jet (40 and 60 m/s) cases, 
the waves shorter than wavenumber 5 are primarily baroclinic 
in nature, while the very long waves have combined baro- 
clinic and barotropic properties. For the stronger (80 m/s) 
jet, the wave becomes more barotropic at all wavenumbers. 
In the real data case, the growth rate curve was approxi- 
mately 4 per day, with similar tendencies to that of the 
NORAPS cases. One possible explanation could be the fact 
that the real data case was the most barotropically unstable 
in the horizontal, and the NORAPS forecast case the least. 
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As discussed in Chapter III, the size of the cyclone 
(low III in Fig. 3.5) was approximately 1000 km at 290000. 
In less than 6 h between 290000 and 290600, the cyclone 
doubled in size as well as intensity. A perturbation with a 
wavelength of 2000 km at the latitude of the cyclone has a 
wavenumber of about 15. These observations are therefore 
consistent with the growth rate associated with the most 
unstable wavenumber. 

The phase speed (Fig. 5.7) for wavenumbers 7 - 30 
for the NORAPS forecast case closely approximate the speed 
of the cyclone. The phase speeds depicted are those for the 
maximum surface pressure perturbation. In general, all the 
phase speeds for the three initial conditions are higher 
than what was observed. This could be related to the fact 
that the model has no friction, which would have retarded 
the speed of movement of the perturbation. The phase speeds 


also decrease with increasing wavenumber (with the exception 


of wavenumber 35), consistent with the results of Gall 
(1976c). 
2. Structure 


The NORAPS forecast case will be used to examine the 
structure of the perturbation more closely because the 
cross-section in this case is oriented through the middle of 
the cyclone (Fig. 5.3). In the NORAPS analysis case, the 
cyclone was northeast of the cross-section (Fig. 5.4). As 
discussed before, the cyclone size corresponds approximately 
to a zonal wavenumber 15, which has the maximum growth rate 
for the NORAPS forecast case. This growth rate as well as 
the phase speed also agree very well with observations. 
Therefore, the structure of the perturbation for wavenumber 
15 in the NORAPS forecast case will be presented, although 
it should be pointed out that the structures for wavenumbers 


10 through 20 are very similar. 
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The variations of the geopotential amplitude pertur- 
bation along the cross-section are depicted in Fig. 5.8. 
Two maxima can be identified, one near the surface between 
Stations AHN and AYS and another at approximately 500 mb 
between BNA and AHN. The positions of these two maxima 
indicate a northward tilt with height. These maxima occur 
in the regions where the Richardson number (Ri) is the 
lowest (not shown) and the vertical and horizontal shears 
are the strongest. 

The maximum temperature perturbation is found near 
AHN, with a similar northward tilt with height (Fig. 5.9). 
This maximum occurs near the surface and coincides with the 
strongest baroclinic zone (Fig. 5.2). The perturbation is 
confined to the lowest layers. 

The maximum perturbation in vertical motion 
(Fig. 5.10) occurs at approximately the 700 mb level, prima- 
rily over AHN and southward. The sign of the vertical baro- 
tropic term (discussed in the next sub-section) calculated 
from the model indicates that this maximum perturbation 
corresponds to an area of rising motion. This area also 
coincides with a region of low Ri (not shown) and strong 
horizontal and vertical wind shears (Fig. 5.2). 

In contrast to pure baroclinic disturbances, these 
temperature and vertical motion perturbations are not coin- 
cident everywhere, with the temperature perturbation to the 
north and at a lower level than the vertical motion pertur- 
bation. These patterns suggest that warm air is rising, and 
cold air is sinking beneath 700 mb (vertical motion values 
increasing with maximum temperature perturbation values). A 
reversal above this level occurs as the temperature pertur- 
bation values decrease with constant vertical motion values. 
This suggests that baroclinic conversion contributes to the 
development of the cyclone at lower levels, but inhibites 


growth at mid-levels. This result is consistent with that 


76 


O HoOdur (1984). These ideas will be discussed further in 
the next sub-section. 

For wavenumbers 10-20 in the NORAPS analysis case, 
the maximum height, temperature and vertical motion pertur- 
bations were located near the anticyclonic side of the jet 
level (coincident with the strong horizontal shear region 
depicted on Fig. 5.1). In these experiments, the barotropic 
mode appears to dominate over the baroclinic mode, resulting 
in the maximum perturbation growth being at higher levels. 
The relative contributions from the barotropic and baro- 
clinic processes will be discussed further in the next 
subsection. At higher wavenumbers (greater than 25), the 
maximum height and temperature perturbations were found near 
the surface, similar to that found for wavenumbers 10-20 in 
the NORAPS forecast case. (Similar results were also found 
Bor the higher wavenumbers in the forecast case). 
Furthermore, in the forecast case, the strongest horizontal 
and vertical shears were located below 800 mb. The strength 
and position of both the horizontal and vertical shears 
appear to dictate the location and amount of perturbation 
growth. Therefore, the emergence of the most unstable baro- 
Elimic or) barcimópa cimode seems | to depend on the relative 
magnitudes of the vertical and horizontal shears. 

The horizontal variations of the geopotential ampli- 
tude perturbations for wavenumber 15 for the NORAPS forecast 
case are depicted in Figs. 5.11 and 5.12. These figures are 
constructed from the profile of the height perturbations in 


Fig. 5.8, using the following equation: 
H(x,y,p)=HO(y,p)sin(2Wx/L + 0) (5.6) 
where L is the wavelength corresponding to a particular 


zonal wavenumber, @ the phase angle and HO(y,p) the maximum 


amplitude of the height perturbation. The vertical profile 
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depicted in Fig. 5.8 is placed at the middle of the plot for 
convenience. Since this is a linear model, the phase angles 
are all relative. 

The maximum surface height perturbation (Fig 5.11) 
is found underneath but somewhat southward of the jet, as 
might be expected. A distinct westward tilt with height of 
the geopotential perturbation can be seen by comparing Figs. 
Ja van 5.123 The height perturbations at 500 mb are 
oriented in a northwest-southeast direction south of the jet 
(Fig. 5.2), which indicates that the horizontal tilt of the 
500 mb heights is opposite to the horizontal wind shear 
south of the jet and parallel to that north of the jet. 


These observations suggest that baroclinic instability 
exists underneath the jet, and barotropic instability 
prevails to the south of the jet. On the other hand, baro- 


tropic stability is suggested to the north of the jet. 

In summary, for the NORAPS forecast case, the hori- 
zontal and vertical structure of the perturbation identified 
by the linear stability model to have the highest growth 
rate bears close resemblance to 'the cyclone. As will be 
discussed further in the next subsection, the barotropic 
modes in the NORAPS analysis case dominated the growth rate 
Spectrum (except for the higher wavenumbers) since the 
maximum perturbations were located near the jet level. 
Although baroclinic instability may also have existed and 
contributed to the growth of the cyclone ina different 
region, this model is not capable of identifying this 
"secondary" instability mode. On the other hand, the 
strongest horizontal and vertical shears in the NORAPS fore- 
cast case were present at lower levels and may therefore 
have allowed the baroclinic mode to become dominant at the 


wavenumbers close to those observed. 
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3. Energetics 
a. Energy Budget Wavenumber Distribution 


Since the zonal flow within the linear model 
consists of a zonal flow and a single zonal wave, the 


kinetic energy equation of the wave can be written as; 








aK ais ou tanp + Eu TA 

ere Uv agr a  (Ú.uv) - uv dp - we (A 
where Ke is the total kinetic energy of the wave, R the 
radius of the earth, Q is the latitude and U is the zonal 
wind. The primed variables represent the eddies. The first 


two terms on the right-hand-side of the equation are 
referred to as the horizontal barotropic conversion term. 
The third term is referred to as vertical barotropic conver- 
Sion. This term is neglected in quasi-geostrophic models, 
and therefore represents a nongeostrophic effect (Gall, 
Oia Jos Hodur (1984) considered this term as a baroclinic 
term because of its relation to the horizontal temperature 
gradient via the thermal wind equation. Orlanski (1968) 
noted that this term may become important for shorter waves 
when the Richardson number is small. Ticetastmeerm repre- 
sents the baroclinic conversion. 

Previous studies by Gall and others have shown 
that waves derive most of their energy by baroclinic 
processes. Shorter ane e greater than 5) are 
barotropically damped in the horizontal. However, the 
vertical barotropic conversion contributes to wave growth 
especially for short waves (wavenumbers greater than 10). 
Therefore, perturbations with higher wavenumbers can extract 


energy from the zonal flow through both baroclinic and baro- 


tropic processes. Hodur (1984) found that for stronger 
jets, the vertical barotropic term increased in importance 
relative to the baroclinic term. On the other hand, Brown 
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(1969) and others found that the horizontal term contributes 
only for the very long waves, and thus have identified the 
planetary waves as being barotropic. 

For the analytic initial conditions with values 
of UMAX set at 40 m/s, the energetics resemble those found 
by Gall ONG With a 60 m/s jet, the results of Hodur 
(1984) were essentially reproduced. Specifically, the baro- 
clinic conversion provides the largest contribution to the 
wave growth, but the vertical barotropic term is more impor- 
tant for higher wavenumbers in the 60 m/s case than in the 
40 m/s case. Except for low wavenumbers, the horizontal 
barotropic term provides a damping effect on wave growth. 

The results for the 80 m/s jet are Shown “In 
Fig.) POR The horizontal term exhibits similar tendencies 
to the experiments with weaker jets, although this term is 
the most dominant at lower wavenumbers. However, the baro- 
clinic and vertical barotropic conversions dominate for the 
higher wavenumbers. For all wavenumbers the vertical baro- 
tropic term contributes to wave development. The baroclinic 
conversion inhibits wave growth for lower wavenumbers but 
enhances shortwave development. With a 6/7 m/s jet, Hodur 
(1984) found the baroclinic and vertical barotropic conver- 
sions to be of equal magnitude. He therefore labeled his 
perturbation a mixed baroclinic wave. 

In the set of experiments using the NORAPS anal- 
ysis as initial conditions, the horizontal barotropic term 
is the most dominant term for all waves with wavenumbers 
less than 25 (Fig. 5.14). The baroclinic term damps wave 
growth for this same range of waves but becomes an important 
contributor beyond wavenumber 25. The vertical barotropic 
term contributes for all waves, especially for higher waven- 
umbers. The large contribution from the horizontal baro- 
tropic term is a result of the barotropically unstable 


initial profiles TES i). 
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It appears that between wavenumbers 25-30, a 
transition from the horizontal barotropic to the vertical 
barotropic and baroclinic modes exists (Fig. 5.14). This is 
consistent with height perturbations calculated from the 
NORAPS analysis case (not shown), where perturbations 
greater than 25 were very similar to those from the NORAPS 
forecast case between wavenumbers 10-20. This transition at 
higher wavenumbers apparently results from the increasingly 
Important contribution from the vertical motion term for 
disturbances of smaller spatial scales. 

For the less barotropically unstable initial 
conditions (the NORAPS forecast case), all terms contribute 
to wave growth, with the vertical barotropic term contrib- 
uting the most to the wave (Fig. 5.15). Furthermore, the 
contribution from the vertical barotropic term increases 
until wavenumber 30 and then decreases. However, the baro- 
clinic term increases andthe horizontal barotropic term 
decreases at higher wavenumbers. 

Therefore, it appears that the structure of the 


perturbation is dependent upon which mode is dominant at a 


particular wavenumber. At those wavenumbers where the 
vertical barotropic and baroclinic terms are dominant, the 
perturbation is confined to lower levels. ` On the other 


hand, the wave will be evident at higher levels if the hori- 
zontal barotropic mode is the strongest. For example, in 
the NORAPS analysis case, the baroclinic mode emerges beyond 
wavenumber 25 when both the vertical barotropic and baro- 
clinic terms become more important. The baroclinic mode 
probably also exists at lower wavenumbers. However, the 
model is unable to identify this mode because of the much 
larger instability from the horizontal barotropic mode. 
b. Vertical Distribution 
The NORAPS 24-h forecast case will be used to 


examine the vertical distribution of the kinetic energy 
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budget terms since the structure of the most unstable mode 
resembles the closest to that observed. In contrast to Gall 
(1976a,c), the eddy kinetic energy (Fig. 5.16) is found to 
be a maximum at 850 mb for wavenumbers 10, 15, 20 and 30 
with a slow increase at the surface for successively higher 
wavenumbers. The fact that the maximum is not found at the 
surface could be attributed to the friction inherent with 
the initial profile. (The maximum -occur at the surface in 
all analytic cases where no friction was present.) Note 
that the maximum eddy kinetic energy is found at the level 
of maximum vertical wind shear. 

The vertical distribution of the eddy kinetic 
energy budget terms for wavenumber 15 in the NORAPS 24-h 
forecast case is depicted in Fig. 5.1/7. The vertical baro- 
tropic term is the largest contributor to wave growth. The 
maximum contributions from the baroclinic and vertical baro- 
tropic terms are largest at 850 mb, whereas the horizontal 
term contributes most at 700 mb. On the other hand, baro- 
clinic damping occurs below 900 mb and above 600 mb, consis- 
tent with the temperature and vertical motion perturbations 
discussed in the previous section. 

Thus, the perturbation can be largely explained 
by the vertical barotropic process at lower levels, with 
baroclinic instability playing a less significant role. The 
baroclinic term actually inhibits growth above the mid- 
levels and at the lowest levels. This finding is similar 
with Hodur (1984) who noted that baroclinic damping occurred 
at mid-levels for waves with wavenumber 40. The shallow 
nature of these waves is also similar with the vertical 


extent of the cyclone in this case study. 
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VI. CONCLUSIONS 


Previous Studies have suggested a variety of mechanisms 
and processes to explain the rapid development of some 
extratropical cyclones. It is reasonable that some favor- 
able interaction between these different processes could 
accelerate the developmental process. This study has iden- 
tified two of these processes as being important for one 
case of rapid cyclone development. 

A detailed synoptic evaluation of the rapid development 
of the Carolinas storm of March, 1985 suggests that the 
atmosphere was probably both baroclinically and barotropi- 
cally unstable during the genesis stage. Therefore, it is 
hypothesized that a combination of baroclinic and barotropic 
instabilities could explain the behavior, structure and 
growth rate.of this cyclone. 

A linear stability model was used to test this 
hypothesis. Three different sets of initial conditions 
along a vertical cross-section perpendicular to the jet were 
employed. These include a temperature analysis using rawin- 
sondes, a NORAPS temperature analysis, and a NORAPS 24-h 
temperature forecast. Geostrophic winds were calculated for 
each of these three temperature fields. The growth rates 
and phase speeds of the most unstable modes were consistent 
with observations. However, the structure of the perturba- 
tion depends on whether the most unstable mode is primarily 
baroclinic somieeba Hors cols The horizontal barotropic mode 
was most evident at upper levels, whereas the baroclinic and 
vertical barotropic modes were found near the surface. The 
unstable mode that best explains the Carolinas storm appears 
to be a mixed baroclinic/vertical barotropic mode. 

This case study of rapid cyclone growth and development 


has shown that the general behavior, Structure and 
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entrseeties of the cyclone is consistent with the results 
from a linear stability analysis. Thus, it appears that 
linear theory can be used to explain the genesis and initial 
growth of the Carolinas storm. Furthermore, baroclinic 
instability alone cannot explain the growth and development 
of the cyclone. In fact, the growth of the cyclone has been 
found to be the result of barotropic and baroclinic 
instabilities. These instabilities can account for the 
initial. growth and development of the cyclone without 
including the effects from convection. 

Research should be continued to further our under- 
standing and prediction capabilities of rapid cyclone devel- 
opment. Specifically, our understanding of barotropic 
instability from both theoretical and applications stand- 
points must be improved. How baroclinic and barotropic 
instabilities interact with each other as well as other 
processes need to be better understood. À more coherent 
picture of how convective and boundary layer processes, 
topographic influences and jet streak dynamics interact and 
contribute to rapid cyclone development should also be 
addressed. The variety of physical processes that are 
present in the rapid development of maritime extratropical 
cyclones suggests that the interaction between different 
scales of motion could be very important. For example, 
Thorpe and Emanuel (1985) have proposed that moist symmetric 
instability, which operates on scales of 100 km and has been 
used to explain frontal structures, can contribute to 
cyclone growth from potential vorticity arguments. These 
ideas need to be explored through observational studies 
uSing a dense data network that can spatially and temporally 
resolve motions on the mesoscale. Numerical and theoretical 
Studies are also necessary to provide further physical 


insights into the problem. 
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